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‘We  have  begun  an  analysis  of  the  optimal  design  of  structure 
in  mechanics  and  biomechanics,  with  emphasis  on  materials  which 
do  not  satisfy  Hooke’s  law.  Our  research  has  been  partly  numerical 
and  partly  analytical;  the  optimal  design  problem  is  highly  nonlim  ar 
with  respect  to  variations  in  the  geometry,  and  this  has  meant  tha 
the  construction  of  efficient  algorithms  is  unusually  difficult. 

They  depend  on  the  formulation  of  the  optimization  problem,  anoT/we 
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20.  continued 


have  made  progress  in  this  area — to  replace  a  nonconvex  cost  fiinction 
by  a  convex  cost.  In  the  case  of  multiple  loadings  our  new  formulation 
is  polyconvex,  and  may  be  the  first  case  of  its  type  to  yield  explicit 
minimization  of  functionals  of  a  Jacobian  matrix. 
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Summary  of  Results 


This  research  has  yielded  progress  in  the  area  of  optimal 
design,  vihlch  Involves  new  mathematical  questions  In  the  formu¬ 
lation  of  a  correctly  posed  optimization  problem.  The  Initial 
statement  of  the  problem  Is  difficult  both  analytically  and 
numerically;  In  many  cases  the  cost  fimctlon  Is  nonconvex. 

Working  Jointly  with  Robert  Kohn  we  have  found  an  equivalent  convex 
problem  for  the  case  of  a  single  loading,  and  our  current  research 
has  led  to  the  calculation  of  an  equivalent  polyconvex  problem 
for  the  case  of  multiple  loadings. 

Our  results  In  these  different  areas  will  be  summarized 
separately: 

1)  For  cylinders  in  torsion  or  in  shear,  we  have  determined 
the  minimal  structure  within  a  given  region  that  can  achieve  a 
specified  torsional  rigidity  or  support  a  prescribed  load.  By 
constructing  explicit  solutions  we  established  a  class  of  problems 
for  numerical  analysis,  and  we  developed  analytical  methods  which 
extend  to  other  geometries  and  other  materials. 

2)  We  prepared  and  tested  a  code  to  link  the  quasi- Newton 
methods  In  nonlinear  optimization  to  the  existing  programs  for 
finite  element  analysis.  Our  package  Includes  two  update  techni¬ 
ques,  the  BFGS  rank  two  update  for  positive  definite  symmetric 
systems  (convex  optimization)  and  the  Broyden  rank  one  update  for 
xinsymmetrlc  systems  with  convective  terms.  Both  have  been  requested 
by  engineers.  Dr.  Engelman  continues  to  work  on  the  development  of 
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a  general  fluid  dynamics-finite  element  code«  which  Is  an 
essential  tool  In  scientific  computing. 

3)  The  max  flow-mln  cut  theorem  of  operations  research  was 
extended  to  the  continuous  case,  with  flows  through  a  domain  and 
local  capacity  constraints  on  the  velocity  field.  A  parallel 
development  by  Irl  has  led  to  applications  In  the  analysis  of 
traffic  flow. 

4)  In  Joint  work  with  Iserles  we  completed  a  study  begun 
much  earlier^  by  establishing  the  upper  limits  on  the  accuracy  of 

a  stable  method  for  hyperbolic  equations.  The  accuracy  Is  restricted 
by  the  choice  of  meshpolnts,  and  high  accuracy  cannot  be  achieved 
In  the  full  region  permitted  by  the  Courant  condition  on  flow 
speeds.  We  studied  both  explicit  methods  and  a  new  class  of 
loQ)llclt  Fade  schemes. 

5)  Ballistic  penetration  problems  lead  to  regions  of  extremely 
rapid  straining  and  to  loss  of  monotonlcty  In  the  stress-strain 
curve.  In  a  report  at  Aberdeen  Proving  Grounds  we  discussed 
progress  on  the  understanding  of  shear  bonds.  We  reported  also 

on  numerical  experiments  with  the  equation  u^.  »  giving 

a  bounded  solution  with  oscillatory  derivatives  In  the  region 
where  o*  <  0. 

We  separate  for  special  mention  two  continuing  research 
projects  which  have  come  from  the  original  results  on  optimal 
design: 


6)  For  two  or  more  loads,  the  mathematics  of  relaxation 
(or  homogenization)  Is  much  more  difficult:  the  relaxed 
form  Is  no  longer  convex  In  the  stresses  ct, t.  Nevertheless 
it  Is  lower  semlcontlnuous  and  a  minimum  Is  achieved. 

This  problem  was  studied  abstractly  by  Morrey,  who  found  a 
property  that  replaces  convexity  in  the  relaxed  integrand. 

We  have  computed  the  first  exeimple  that  fully  illustrates 
this  theory.  It  arises  in  elastic  design  for  given  compliance 

inf  J  ^^{|a|+(T|?^03  dx  =  min  J*  G(o,t)  dx. 

It  is  the  integrand  G  which  we  fo’jnd,  and  it  is  polyconvex — 
convex  function  of  a,  t,  and  D  =  We  now  have  two  ways 

to  compute  G.  One  uses  the  bounds  found  by  Tartar  on  the 
moduli  of  a  composite  elastic  material,  given  the  fraction  p 
of  the  original  material  used  to  construct  the  composite.  Our 
method  is  to  consider  all  possible  composites  near  each  point, 
so  that 

1  +  |a|^  +  |t(^  -»•  p  +  (Aa,a)  +  (At,t) 

where  A  gives  the  compliance  of  the  composite,  "nien  we 
minimize  over  all  A(p)  allowed  by  Tartar,  and  over  all  p. 

Our  second  method  is  more  genera.'.,  but  the  computations 
are  again  not  easy.  It  finds  the  "polyconvex  hull"  of  the 
graph  of  the  original  integrand. 


We  have  Just  begun  on  these  ideas.  They  are  certain  to 
lead  further— to  a  duality  theory^  If  possible^  and  to  a  deeper 
understanding  of  nonconvex  variational  problems  In 

7)  Jointly  with  the  orthopaedic  research  group  at  Beth 
Israel  hospital*  we  are  studying  Wolff’s  law  for  the  structure 
of  bone  fibers:  their  density  and  direction  are  determined  by 
the  imposed  loads .  The  first  niimerlcal  experiments  show  agreement 
between  the  features  of  a  theoretically  optimal  design  euid  the 
actual  distribution  within  bones.  We  plan  to  compare  their 
stereological  description  and  the  output  from  finite  element 
analyses  with  the  optimal  designs  computed  from  further  numerical 
optimizations.  The  code  approximates  the  Michell  theory  of 
optimal  trusses,  and  leads  in  two  dimensions  to  orthogonal  nets 
of  Hencky-Prandtl  type.  The  orthogonality  of  bone  fibers  can  be 
studied  experimentally;  if  confirmed,  the  three-dimensional 
analysis  will  lead  to  new  questions  in  optimal  design.  We  believe 
there  are  applications  of  both  6)  and  7)  to  the  construction  of 
composite  materials. 
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